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Abstract 

In yeast {Saccharomyces cerevisiae) and plant roots [Arabidopsis thaliana) zinc enters the cells via influx 
transporters of the ZIP family. Since zinc is both essential for cell function and toxic at high con- 
centrations, tight regulation is essential for cell viability. We provide new insight into the underlying 
mechanisms, starting from a general model based on ordinary differential equations and adapting it to 
the specific cases of yeast and plant root cells. In yeast, zinc is transported by the transporters ZRTl and 
ZRT2, which are both regulated by the zinc-responsive transcription factor ZAPl. Using biological data, 
parameters were estimated and analyzed, confirming the different affinities of ZRTl and ZRT2 reported 
in the literature. Furthermore, our model suggests that the positive feedback in ZAPl production has 
a stabilizing function at high influx rates. In plant roots, various ZIP transporters arc involved in zinc 
uptake. Their regulation is largely unknown, but bZIP transcription factors are thought to be involved. 
We set up three putative models: activator only, activator with dimcrization and activator/inhibitor. 
These were fitted to measurements and analyzed. Simulations show that the activator/inhibitor model 
outperforms the other two in providing robust and stable homeostasis at reasonable parameter ranges. 

Author Summary 

The heavy metal zinc is both an essential micronutrient for all living organisms and a toxin at high 
concentrations. Zinc uptake in plants has become an especially interesting research area, because certain 
hyperaccumulating plant species can be used to clean up zinc-contaminated soils in mining or industrial 
areas. On the other hand, crops may be manipulated to accumulate zinc in order to overcome nutritional 
zinc deficiency, which is a problem in many developing countries. Since plants and fungi rely on the zinc 
content of the surrounding medium, they need to have mechanisms providing a tight uptake regulation at 
varying external zinc concentrations. In view of the promising applications it is important to understand 
these mechanisms. Here, we want to support experimental findings with some purely theoretical insights. 
We use quantitative data from plants and yeast to develop feasible models, analyze their properties in 
a general framework and discuss their biological relevance. The results confirm experimental hypotheses 
and may give input for further measurements. 



Introduction 

Zinc is a heavy metal and micronutrient that plays an important role in all living organisms and is 
particularly essential for the growth of higher green plants [1]. It is part of the functional subunits 
or cofactor of more than 300 proteins, among them the class of zinc-finger-proteins as well as RNA- 
polymerases. In addition, it has been reported to protect plant cells from oxidative stress mediated by 
reactive oxygen species (ROS) [2] and may act as an intracellular second messenger [3]. 

In higher doses, however, zinc becomes toxic. Toxicity is far less frequent than deficiency, but likely 
in plants growing on contaminated soils, e.g. in mining or industrial areas. Most plants react to elevated 
zinc levels with toxicity syndromes, such as reduced growth and leaf chlorosis [4]. Only specialized 
zinc-hyperaccumulating species are able to tolerate high levels without impairment [5]. In order to do 
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so, they possess mechanisms for both the increased uptake of zinc from the soil and its sequestration 
and detoxification [B]. These mechanisms are subject of ongoing research, as they implicate interesting 
applications in phytoremediation or nutritional enhancement [7] . 

Avoiding both deficiency and toxicity, plants need to take up their required amounts of zinc. Unlike 
animals they cannot adapt their nutrition accordingly, but depend on the zinc content of the soil. This 
content may vary considerably in different locations and under different conditions. How are plants able 
to adapt to this variety? 

Charged zinc ions are unable to cross cell membranes freely [8]. Instead, they are taken up by 
specialized transporter proteins. To provide a sufficient intracellular zinc concentration without reaching 
toxicity, these transporters need to be tightly regulated. The regulatory mechanism has to consist of two 
parts: sensing of the intracellular zinc concentration and reaction to changes by controlling the amounts 
of zinc transporters. 

The sensing of changes in zinc concentrations must be very sensitive, because the actual internal zinc 
concentration is believed to be very small. Zinc ions bind to various intracellular proteins, are chelated 
and sequestered into specific cellular compartments, such as the vacuole [9]. Therefore, although the 
zinc content in the entire cells may be in a millimolar range, the actual concentration of free zinc ions is 
estimated to be in a femtomolar range [10| . Zinc influx carriers are thought to be regulated by this pool 
of free zinc ions plus ions that are loosely bound to chelator proteins and can be set free to bind to other 
proteins with higher affinity. 

Models of homeostasis 

Homeostatic regulation in biological systems is based on genetic regulatory systems, and ultimately, on 
concentrations. These are positive, which constrains the possibilities of control substantially. [H] showed 
that the positiveness constraint of a perfectly / robustly regulating enzyme leads to the need of two 
separate control mechanisms: for influx and efflux. The homeostatic model proposed by is 



where S is the regulated species, / = I{S, R) > and E = E{S, R) > are the influx and efflux, respec- 
tively, R is the regulator, fc is a coefficient (not necessarily positive) and Ss is the set point concentration. 
The above model results in non-physical negative concentrations of the regulator Independently 
of the type of mechanism sought after, the negative term in dR/dt needs certain properties to achieve 
robustness based on positive concentrations. The approach is to have a term which is linear in R for 
small R (positiveness), but becomes almost independent of R for larger R (robustness) [TT] . 

Eq. ([1]) is an oversimplification of homeostatic control in cells, as substantially more complex mecha- 
nisms are needed (compare Figs. [T]and|4]). Also the concept of perfect control is an idealization. Control 
of zinc fails in cells for low and high external conditions. The presence of oscillations in perfect homeosta- 
sis, [12], poses a problem to living organisms. Strong osciflations could lead to transient, very high and 
potentially lethal concentrations. Prescinding from perfect regulation could be a compromise between 
avoiding strong bursts and achieving good control. 

Based on biological information available, we will develop several putative models of influx homeosta- 
sis in plant root cells. In the first part, a general influx regulation model based on an ordinary differential 
equation system describing transporter gene expression, will be developed and non-dimensionalized. Us- 
ing the general model, the biological model for yeast of [13] will be translated into a corresponding 
mathematical model. This model is simplifled and fitted to transcript level data via a non-linear opti- 
mization method |14j . The mathematical properties of the steady state, such as stability, are analyzed 
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and discussed. In the second part, the experiences won with the yeast model are used to pose three 
models for plant roots. The possibilities are manifold, for which reason we restrict the models to the 
most simple cases of: activator only, activator with dimerization and activator-inhibitor. 



Results and Discussion 
General model 

The zinc homeostasis mechanisms presented in this manuscript can be arranged into a general model, 
which will be developed in this section. Zinc homeostasis can be split into two components: short and 
long term regulation. Short term regulation is fast but rough, while fine tuning is done by long term 
regulation. The time scale of short term regulation is less than two hours in plant roots |15| . Long term 
regulation has a substantially larger time scale of several hours, days, weeks, etc. 



We are interested here in short term regulation, which is local in the sense that the processes occur 
at the level of single cells in plant roots. Signals besides the fluxes seem not to be transmitted between 
cells or tissues. This is of course not the case for long term homeostatic control, which will rely definitely 
on signals transmitted from tissue to tissue. The short term response in plant roots and yeast cells is 
expected to follow similar laws and can be subdivided into the phases 
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The zinc status is measured in the sensing phase, decisions are taken in the transduction phase and 
changes in cytosolic concentration occur in the reaction phase. As mentioned in Sectiouo both infiux and 
efflux can be adapted to achieve homeostatic control. In plant roots as well as in yeast cells, adaptation 
of the expression of influx transporters poses the major component of zinc regulation |151ll6j . 

Based on the concept presented in Eq. 1^, the models considered in this manuscript have the following 
structure 
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Reaction: 
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where Z and Z"^ are the cytosolic and external zinc concentrations, respectively, Ai are activators, li 
inhibitors, Ti and Ei influx and efflux transporters, respectively, Gi and Mi the levels of gene expression 
and mRNA of T^, respectively, and pAi and pji arc model dependent production terms. The total 
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activation and repression are 



J = l j,k=l J = l 



and the function f{Z, K) describes saturation of the transporters 
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Sensing is assumed to take place at the level of the activators Ai and inhibitors /; . The possibility 
that the transporters Ti sense the cytosolic zinc concentration Z directly was also introduced. To achieve 
regulation, the total activation Ai has to decrease with higher Z values (see Section^. Transduction is 
modeled in the usual way [17]. Three equations per protein are needed, namely for: gene activity Gi, 
transcription into M, and translation into Ti. The activators are introduced as essential transcription fac- 
tors, which activate the gene transcription. The quadratic form in Eq. (j4]) allows to include dimerization. 
The inhibitors inhibit either the activators or repress through li directly gene activity. Gene repression 
was assumed to be non-competitive and fast compared to activation, i.e. it is in quasi-equilibrium and 
Kij are equilibrium constants. The reaction phase is described by an equation for the cytosolic zinc 
concentration, which contains essentially the difference between influx and efflux mediated by Ti and Ei , 
respectively, and a transporter independent consumption -7Z. Regulation of the efflux transporters Ei 
was left out of Eq. ([3|), as these vary only slightly in roots and no information on yeast was available. In 
essence, these proteins would follow a similar transduction system as the influx transporters Ti. 

Non-dimensionalization of transduction in Eq. ([3]) is straightforward using 
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Non-dimensionalization of the sensing equations depends on the particular structure of the production 
terms. The decay terms can be non-dimensionalized choosing 
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Yeast 

The regulation of zinc uptake in yeast cells (Saccharomyces cerevisiae) has been studied in much detail and 
found to be a combination of two systems with high and low affinity for zinc ions. A similar distribution of 
high and low affinity transporters has also been found in wheat plants [18] and is thought to exist in other 
plants as well }19j . A schematic overview of the system can be seen in Fig. [TJ Zinc ions arc transported 
with high affinity by ZRTl (zinc-responsive transporter) and with low affinity by ZRT2, which both 
belong to the ZIP (zinc-, iron-permease) family. ZRTl has been found to be strongly regulated by the 
intracellular zinc concentration and almost exclusively active under conditions of zinc deficiency |20) . 
ZRT2 has been reported to guarantee a basic zinc uptake level under normal zinc-repletc conditions |21| 
while being repressed under zinc deficiency [22| . Further studies have shown that both ZRTl and ZRT2 
are activated by the transcription factor ZAPl (zinc-dependent activator protein) jl3j . which binds to so- 
called zinc responsive elements (ZREs) in the promoter regions of the respective genes. Under conditions 
of elevated zinc concentrations, the activity of ZAPl is reduced and production of ZRTl and ZRT2 
decreases. Inactivation of ZAPl occurs most likely by direct binding of free zinc ions, although further 
signaling molecules may also be involved in this process. By binding to its own promoter region, ZAPl 
regulates its transcription introducing a positive feedback mechanism and presumably allowing an even 
stronger response to zinc- limiting conditions, |16| . In addition to the transcriptional regulation, ZRTl 
is also regulated by a post-translational mechanism |16] . While it is a stable membrane protein under 
zinc deficient conditions, ZRTl is ubiquinated and subjected to endocytosis for high intracellular zinc 
levels. The exact details of this mechanism have been investigated by [23], but it is yet unknown whether 
zinc ions bind directly to ZRTl to induce its ubiquitination, or whether other zinc-binding proteins are 
involved. It has been proposed that the combination of transcriptional and post-translational regulation 
allows for a very quick response to changing environmental conditions and thus prevents a toxic zinc 
shock, [16j . An overview of the yeast zinc uptake mechanism is presented in Fig. [TJ 

Model 

As described above, zinc uptake regulation in yeast comprises the two zinc transporters ZRTl and ZRT2, 
as well as the transcription factor ZAPl as the only activator, which is directly inhibited by zinc ions 
without an inhibitor. The production of the activator, which corresponds to the term pAi{Ai, ...) in the 
general model Eq. ([3]), is a system of sensing, transduction and regulation by itself, because ZAPl acts as 
its own transcription factor through a positive feedback loop. While ZRTl is simply activated by ZAPl, 
ZRT2 is both activated and repressed by the same molecule [22]. Therefore, we assume a model with two 
binding sites of ZAPl close to the ZRT2 gene, one activating and one repressing. The total inactivation 
li (see Eq. ([5])) introduces this mechanism into the general model Eq. ([3]). Here, the inhibitor is equal to 
the activator and only the ZRT2 gene is affected: Zi = and I2 = K'2 A. 

Following the syntax of the general model and the non-dimensionalization derived in Section „ we 



Modeling zinc uptake regulation 



6 



obtain the following system 
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The post-translational regulation of ZRTl is given by the term — 7n Tti Ti Z. For simplicity the term 
—7 Z accounts for all zinc consumption processes. These may include export from the cell through 
zinc efflux transporters, sequestration into the vacuole and other compartments as well as binding and 
chelation of zinc by various proteins in the cytoplasm. 

The trivial solution (all species zero) is a steady state of Eq. (|6]). There is at least one non-trivial steady 
state, which for the activator ZAPl can be written as a function of the intracellular zinc concentration 

A^^^-' (7) 
1 + TaZ Ka 

For A not to become negative, this equation poses the condition Ka > 1 + F^ Z, which implies that for 
large Z the non-trivial and trivial solutions cross. A detailed analysis of this case is presented below. 
The case of total deficiency brings insight into some of the parameters. We find as expected Z ^ for 
Z'^ — > 0, which means that A — 1 — 1/Ka- From biological point of view, A is expected to shoot to a 
value near to 1 for total deficiency, which has as a consequence that Ka ^ 1. Assuming that A « 1 for 
Z'^ 0, the concentrations of the transporters Ti and T2 behave as 

J\ —, and ■ ^—r-, 7 for Z" . 

l + l/Ki l + llK2 + K'^/K2 + K'^ 

High affinity of ZRTl and low affinity of ZRT2, i.e Ti w 1 and T2 w for ^ 0, are obtained when 
the conditions A'l > 1 and K'^ + K'^jK^ + ^jK^ > 1 are fulfilled. Considering K2 ~ K^ \, the 
second condition is essentially K2 ^ 1. Expression of ZRT2 is maximal for a ZAPl concentration of 
A ~ (i^2 ATj)^^/^, while expression of ZRTl rises monotonically with A and reaches its highest value 
for Z'^ — ^ 0. For a given activation 1^2 , repression K'2 has to be large to shift the expression maximum 
towards low A and high Z'^ . 

Using the quantitative data measured by |13j and j22j . we estimated the model parameters via opti- 
mization. These measurements are stationary, and thus, the system could be reduced into one with the 
four unknowns A, Ti, T2 and Z. The parameters obtained are listed in Table [1] These reflect perfectly 
the above conditions for KA^ Ki, K2 and ATj- The model reproduces very well the measurements (Fig. 
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Roles of ZRTl and ZRT2 

[13] proposed that ZRTl and ZRT2 play different roles in zinc uptake of yeast cells. While ZRTl is 
most active only in zinc-deficient cells, ZRT2 is transiently active also in zinc-replete cells with external 
zinc concentration around 1000 fiM . This implies that under low external zinc concentrations ZRTl 
contributes the most to the overall zinc uptake, while under high external zinc concentration, ZRT2 acts 
as the major transporter. Such behavior is confirmed by our model. 

Fig. |3] A presents the relative contributions the the total flux. At low external concentrations ZRTl 
is responsible for about 80% of flux, while at replete conditions (above 500 /iM) ZRT2 dominates. ZRTl 
seems indeed to act as a high affinity transporter with a Michaclis constant K{ = 139 ^M, while ZRT2 has 
less affinity reflected by a substantially larger K2 = 2584 fiM. A similar ratio was found by [21], although 
their values are several orders of magnitude lower. This discrepancy comes from the assumption made 
in |21j that the mechanism is based on Michaelis-Menten kinetics and their low values are reproduced 
when Michaelis-Menten is fitted to our simulations. The affinity of the ZRTl and ZRT2 systems are 
not completely determined by K[ and K2, respectively. These constants have to be larger than the 
optimal concentration of the corresponding system, as saturated transporters cannot pass information on 
external zinc status {f{Z^,Kj) 1 — const for ^ K*). The optimal concentration for ZRTl is at 
total deficiency, while ZRT2 is most active at 430 //M (Fig. [2]). 

A strong repression of ZRT2 is essential to achieve a maximal expression at high external zinc concen- 
trations (see Table [ij. However, a strong repression results also in lower gene activities, which explains 
why ZRT2 has a much lower expression level than ZRTl (Fig. [2] and j22]). To counteract the lower 
expression level, ZRT2 needs either to transport zinc at higher rates or more copies need to be produced. 
This is reflected by the coefficient k, which suggests that ZRT2 is six times more effective in transporting 
zinc than ZRTl. Assuming that ZRTl and ZRT2 have similar transport rates, k pa 6 could be a hint for 
posttranslational regulation of ZRTl, although direct posttranslational regulation via Vti was shown to 
be not signiflcant here (F-test: P > 0.05). 

ZRTl and ZRT2 were found to be activated equally well by ZAPl, as reflected by the small insignifi- 
cant difference between Ki and K2- The self-activation constant Ka of ZAPl, is four times smaller than 
Ki and K2- This suggests that ZRTl and ZRT2 have four promoters instead of one in ZAPl, which is 
in concord with the experimental results of |13j . 

ZAPl transcriptional feedback 

The feedback loop generated by ZAPl acting as its own transcription factor brings interesting properties 
into the model. |16| proposed that this feedback allows a stronger reaction to zinc-limiting conditions. 
Our model suggests that the advantage is rather for zinc-replete conditions. The steady state Eq. ([7]) of 
ZAPl becomes negative for Z > (Ka — i)/^A ~ 0.15 and crosses the trivial steady state. Unless these 
two steady states exchange their roles, the model would become non-biological at the bifurcation. Based 
on the fitted parameters, the bifurcation is normally reached at very high external zinc concentrations. 
To examine the behaviour of the model at the bifurcation, we introduced a ZRTl and ZRT2 independent 
path into the cell. The path could for example be another transporter not regulated by ZAPl and 
shifts the bifurcation towards lower Z'^. Without considering any details of these processes, the simplest 
modification is to include an additional constant zinc infiux term az to the last equation in Eq. 
The bifurcation is illustrated in Fig. [3] B. There are at least two steady states, where one is trivial 
(A = Ti = T2 = and Z = az) and the other is positive for small az- The stability of these are 
exchanged at the bifurcation. For low az the positive steady state is stable, while the trivial steady state 
is unstable. When the steady states cross at the bifurcation, the trivial solution becomes stable while the 
now negative steady states becomes unstable. The positive steady state is literally trapped by the trivial 
steady state. From the biological view the ZAPl feedback allows the system to completely switch off 
expression of ZAPl and thus of ZRTl and ZRT2. In a mechanism without feedback, ZAPl expression 
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would just fall asymptotically to zero for increasing zinc influx. Therefore, we conclude that the feedback 
of ZAPl is advantageous for zinc-replete conditions. 

Plant roots 

In plants, zinc is taken up from the soil and transported into the root cells. Unlike in unicellular organ- 
isms, zinc needs to be transported into further tissues: xylem, stem, leaves, etc. A number of different 
transporter proteins are involved. There are three families of transporters for zinc: ZIP, HMA (heavy- 
metal- ATPases) and MTP (metal tolerance protein) or CDF (cation diffusion facilitator). Members of 
the ZIP family are believed to act as influx carriers, including uptake of from the soil (similar to ZRTs in 
yeast). HMAs accomplish efflux of zinc, e.g. from roots into xylem vessels, while MTPs are involved in 
sequestration into compartments, such as the vacuole |23]. The main root influx transporters are ZIPl, 
ZIP2, ZIPS, ZIP9, and IRT3 [25], while ZIP4 locahzes to chloroplast [19]. These transporters are highly 
expressed under conditions of zinc deficiency, whereas their expression decreases quickly when zinc is 
added to the media [TSj. The exact mechanism of this regulation is still unknown. Recent results showed 
that at least ZIP4 in Arabidopsis thaliana is regulated by transcription factors of the basic-region leucine 
zipper (bZIP) family: bZIPlQ and bZIP23, [55]. These factors bind to a ZDRE (zinc deficiency response 
element), which found in the upstream regions of the ZIPl, ZIP 3, ZIP 4, and IRT3. 

Unlike the ZAPl transcription factor in yeast (see Sectioup), bZIP19 and bZIP23 transcription factors 
do not have a zinc binding site [26]. It is unclear how these sense the intracellular zinc status. Existence 
of further players that bind zinc and act as inhibitors of bZIP19 and bZIP23 have been proposed [27] . 
Transcription factors of the bZIP family have been studied in other regulatory networks and are known 
to be regulated post-transcriptionally in various ways [28) . Generally, bZIP transcription factors and in 
particular bZIPlQ and bZIP23 act as dimers [29]. They are partially redundant [26] and it is believed 
that they preferentially form homodimers, but may also interact to constitute heterodimers [30j . 

Our model focuses on the specific transition of zinc into the root cell space. By restricting the model 
to this specific situation, a similar approach as the one for yeast in Section □ can be applied. We start 
with a simple model based on one zinc dependent activator. Hereafter, the advantage of dimerization 
is analysed and a more involved model based on an activator/inhibitor pair is presented. Using the 
measurements of |15j . some of the parameters are obtained via optimization and a F-Test is used to 
conclude which model is the best. Finally, we analyse the relation between stability and robustness of 
the activator/inhibitor model. 

Activator 

Here, we assume that regulation takes place by one zinc dependent transcription factor (see Fig.|4] (i) for 
a scheme). Based on the general model Eq. ([3]), we set ua = tit = 1 and n/ = and avoid unnecessary 
notation by dropping indexes (e.g. A = Ai and K = Kn, etc.). Sensing is assumed to take place only at 
the activator level (/3t = 0). Also the possibility that the activator acts as a dimer is ruled out {a^^j = 0). 
Efflux transporters arc assumed to be non-saturable, allowing to combine efflux/consumption into one 
term —7 Z. In contrast to the case of yeast in Section „ there is no specific information on the production 
of the activator. To keep the system simple, we introduce a constant pool ^0 of activator, which is 
distributed into active and inactive molecules. The net production is set to aA {Aq — A) and 



PA = aA Aq and aA = 7a 



(8) 
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The non-dimensionalized system is 
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Z = jj- (^-i(A' + l)±(A'r^/(ZM^*) + i(A' + lf)*) . 

The steady state with Z < is not biologically relevant and is not considered. For total deficiency 

Z^O and G ^ ^ for Z'' ^ . (11) 

K +1 ^ ^ 

Biology suggests that gene expression will shoot to a very high value and G should be near to one: K ^ 1. 
For replete conditions 



z ^ -i(/v + i) + (A'r^ + i(x + i 



for Z" -^oo , (12) 



where f{Z^,K*) — >■ 1 and A' ^ 1 were used. Biology suggests that gene expression should be small for 
high external zinc concentrations 

Fa > A' > 1 . (13) 

The steady state for varying Z*^ depends on three parameters AT, F^ and AT*. While A'' is a property 
of the transporters, K and Ta determine gene activity for extreme conditions. A value AT* = 13 iiM for 
ZIPl was published by [IS] and used here. Assuming that gene activity reaches at least 95% for total 
zinc deficiency, one obtains 

K >20 . (14) 

Determination of AT from measurements would need data at very low zinc concentrations, which is 
uncertain and was not available to the authors. For this reason, an empirical value of AT = 20 was used. 
Only F^ stays undetermined and a value of 41138 was obtained by fitting the model to published values 
of ZIPS expression [TS] . All parameters are listed in Table [5J 

Fig. [S] shows the steady state for varying Z'^. Gene activity falls slowly for increasing Z"^ resulting in a 
continuously increasing internal zinc concentration. The mechanism breaks for extreme zinc conditions: 
undersupply at low Z*^ and oversupply for large Z"^. The reason for oversupply is the activator reacting 
slowly to changes in Z''. The model offers only a mean to set the extreme gene expressions via AT and 
F^, and does not allow to set the transition steepness between these. F^ is also very large compared to 
the value determined for yeast ('^ 60 times larger; Table [1]), making the model even more unlikely. 
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Dimerization 

The transcription factors bZIP19 and bZIP23 are known to act as dimers [52]. We set a^j ^ and 
aij = in the general model, so that the activator functions only as a dimer. A scheme of the model is 
presented in Fig. |4] (ii). The total activation is here A = K A^, while the rest stays the same as in Eqs. 
© and (|10p . meaning that only gene activity needs to be adapted 

rlC 

— = {K A^l - G) - G) , (15) 

G ^ 



K + (1 + zf 

Because 1 + T a Z < {1 + Ta Z)'^, gene activity is smaller with dimerization than without (compare Eq. 
([T0|) ). Gene activity reacts more sensible to changes of zinc status than in the non-dimersing case (Fig. 
[5]). The transition between gene on and off is steeper, rendering a more robust mechanism. Fitting the 
model to the measurements delivers F^i — 1844, which is ca. 20 times smaller than in the non-dimerizing 
case and substantially nearer to the value for yeast. From an evolutionary point of view, dimerization 
allowed to down-regulate the transporters more strongly with less binding affinity. Also, by assuming 
that the variances of the measured values are proportional to these, one finds that is for the model 
with dimerization less than half as for the one without when fitted to measurements the of [T5]. In total, 
the model with dimerization outperforms the model of Sec. „ although these have the same number of 
degrees of freedom. 

Activator /Inhibitor 

Including dimerization delivered a better fit to the measurements. However, a systematic deviation for 
higher Z'^ was found (Fig. [5]). Following |22l proposition of intermediate steps in sensing, we propose a 
mechanism involving an activator/inhibitor pair. Assume that these interact when they are not bound 
to the DNA, the pairs cannot activate the gene and zinc is sensed only by the inhibitor (Fig. |4] (in)). 
Applying this assumptions to the general model Eq. 1^ delivers ua = nj = ut = 1- As in Sec. „ 
dimerization is included by using the total activation A = K A^ . Production of activator is set as in the 
activator only model (Eq. ([U). Sensing occurs at the level of the inhibitor 

Pi ^ aj Iq Z , ai = 13 j and Pa = . 

Transcription and translation are the same as in the dimerizing activator case. The equation for Z stays 
the same, meaning that the main differences to Eq. ([9|) are 

^ = ja{KA^l-G)-G) , 

A A 

= 7A(l-rA/-^) , (16) 

^ = 7/(r/^-r'A/-(i + r,z)/) . 

If Z is considered to be a parameter in the above system, the steady state is 

G = 

A = 



A' + (i + r/)2 



r7 ' 



\ ( z-c l^_^M ^ 
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where C = T'/TTj and ^ = l/F/. The solution with / < is biologically irrelevant. For totally deficient 
conditions 

7^0, A^l and G ^ — ^ for ^ . 

K + 1 

The case of very high external zinc needs to include the expression for Z. Instead of determining what 
happens for — > c», we determine the behavior for large internal concentrations 

I->i, A^YTf """"^ A- + (f+r)2 

The same biological conditions as those listed in Eqs. and ([Ti]) arc found here. In contrast to the 
activator models, gene activity does not go to zero for Z oo. Again, the constants F and K determine 
gene activity for extreme zinc statuses. The steady state values depend on two more constants: <^ and 
^. The first term in / is zero for Z ^ (. For Z < ( the term is negative and has to be compensated by 
the slightly larger square root term. Is Z > C,, then I{Z) grows fast, inhibits the activator and leads to a 
strong reduction of gene activity (compare Fig.[6]j4). Thus, C determines the internal zinc concentration 
for switching the gene from on to off. The constant ^ determines the steepness of the transition between 
the on and off states (Fig.^B). A small ^ corresponds to a strong binding affinity F/ between zinc and 
inhibitor. The switching steepness is also affected by F, as it weights the first term under the root. Large 
F result in steeper switches (effect similar to decreasing ^; Fig. [7] 5). 

The activator/inhibitor model renders a better and more robust homeostatic control mechanism than 
the activator only models (Fig. [5]). The reason is the steep genetic switch obtained by the inclusion of an 
inhibitor, which reacts strongly to the internal zinc status (Fig.lH]). Fitting the model to the measurements 
delivered F 38 and C, — 4:A ■ 10~^ (Table [2|). ^ cannot be determined by a fit, because a robust 
mechanism is sought after and in that regime the model becomes almost independent of ^ (compare Fig. 
[S]B). Therefore, a value of the same order as Ta for yeast was used {£_ = 10~^ ^ F/ = 1000 and Ta = 714 
for yeast). The model describes the measurements very well (Fig. [S]), which is also a consequence of the 
small number of degrees of freedom. No systematic deviation for large Z*^ was found for this model. An 
F-Test showed that the activator/inhibitor model performs better irrespective of having one parameter 
more (P < 0.05). 



Robustness and instability 

[12] showed that a perfect homeostatic control can lead to undamped oscillations. In the case of a toxic 
compound, oscillations could imply lethal peaks. Therefore, the stability of the activator/inhibitor model 
was analysed. Dynamics and stability depend on the time scales involved in the mechanism. The authors 
could not find published values for these. Similar values to those listed in [31] were used, where the 
products were assumed to decay four times slower than gene activity. The reader should keep in mind 
that the specific choice of the time scales influences stability, but the relation between robustness and 
instability found below should keep its validity. 

F sets the robustness of the activator/inhibitor mechanism by increasing the steepness of the genetic 
switch. This leads also to instability of the steady state and to undamped oscillations (Fig. [7]). The 
oscillation amplitudes for 10 F are shown in the graphs. The model is stable for the nominal parameters 
(Tabled]). Increasing robustness via F lead to instability (Fig. [7]yl). During one oscillation period the 
internal zinc concentration reached up to 3.5 times the steady state value, meaning that strong and 
possible toxic periodic peaks of zinc are produced. These peaks exceed the steady state values of the less 
robust mechanism (Fig.[7]i?). Toxicity for high external zinc conditions could be either because of stable 
high internal zinc concentrations (non-robust mechanism) or due to toxic high amplitude oscillations 
(robust mechanism) . Reducing robustness could be a strategy to avoid strong zinc bursts and cell might 
use other mechanisms to damp strong oscillations, such as buffering and sequestration. 
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Methods 

The ordinary differential equation systems were simulated with either an explicit eighth-order Runge- 
Kutta method or an implicit Rosenbrock stepper for stiff differential equations. Steady states were 
calculated by Newton's method in combination with a path following method for varying parameters. 
Jacobians were calculated analytically. The model parameters were determined by fitting the model to 
measurements. For this purpose, Brent's algorithm was applied to minimize |14il32j . The standard 
deviation of a measurement was assumed to be proportional to its value and the relative error (17%) was 
chosen such to obtain a reduced of ttis order of one. This way, low and high values had the same 
weights and were fitted equally well. Penalties were added to to avoid negative parameter values. 
The confidence intervals were obtained by calculation of the covariance matrix via the Hessian of |32j . 
The measurements in [13l[22] were combined and scaled correctly. Determination of part of the scaling 
factors were included into the fitting process, while the rest was prescribed with given values (personal 
communication of D. Eide). 
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Figure 1. Yeast: scheme of zinc influx regulation model. ZAPl is inactivated by zinc and 
activates transcription of the transporters ZRTl and ZRT2. 



Modeling zinc uptake regulation 



16 




Figure 2. Yeast simulations: Comparison between measurements and simulated steady states of 
ZAPl, internal zinc, ZRTl and ZRT2 for varying external zinc concentration. Measurements: ZRTl 
and ZRT2 by [22], ZAPl and zinc by [13]. 
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Figure 3. Yeast: Role of ZRTl and ZRT2 and ZAPl feedback. A, contributions of ZRTl or 
ZRT2 to the total zinc influx for varying external zinc concentration. B, ZAP activity for varying values 
of ZRT independent influx az- The stable solution is marked with a solid line, the unstable solution is 
dotted. 
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Figure 4. Plant roots: Scheme of the three models of zinc uptake regulation, (i) Activato: 
only, (ii) Activator with dimerization, (iii) Activator/Inhibitor model. 
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Figure 5. Plant roots: Steady states of the different regulation models. The activator only, 
dimerizing activator and activator/inhibitor pair with dimerization. Measurements of [l5] arc also 
shown. 
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Figure 6. Plant roots: Activator/Inhibitor model with dimerization. A, steady state values of 
inhibitor, activator and gene activity in dependence of internal zinc concentration. B, steady state gene 
activity in dependence of internal zinc status for varying ^. Dashed curve corresponds to the nominal 

e = 10-3. 
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Figure 7. Plant roots: Robustness and stability. Robustness and stability of the 
activator/inhibitor model for O.IF, T and lOF and varying external zinc concentration. A, real part of 
largest eigenvalue. B, internal zinc concentration. Minimal and maximal values of limit cycle shown for 
unstable steady state (lOF). 
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Tables 

Table 1. Yeast: parameters 



Parameter Value ± s.d. 



Ka 


109 


±38 


Ki 


450 


±307 


K2 


444 


±119 


K', 


2171 


±1191 


Ta 


714 


±600 


Fti 


29.6 


±31.5 


K 


6.3 


±3.0 


K{ 1 [lM 


139 


±65 


2584 


±1511 



Parameters values and standard deviations obtained 
by fitting the model to measurements published 
in [13] and [22|. 
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Table 2. Plant roots: parameters 



Parameter Act 


■ only 


Act. dimer. 


Act./Inhib. 
dimer. 




13 


13 


13 


K 


20 


20 


20 


Ta 


41138 


1844 




r 






38 


r' 






167.2 








1000 








4.4- 10-3 


e 






10-3 


* Value for ZIPl, 


ED; 


tc = r7rrj; u 


= i/r/. 



Plant roots: parameters used in the simulation of the 
activator only, dimerized activator and the dimerized 
activator /inhibitor models. 



